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Abstract 

We study the behaviour of charge oscillations in Superconducting Cooper Pair Boxes 
weakly interacting with an environment. We found that, due to the noise and dissi- 
pation induced by the environment, the stability properties of these nanodevices differ 
according to whether the charge oscillations are interpreted as an effect of macroscopic 
quantum coherence, or semiclassically in terms of the Gross-Pitaevskii equation. More 
specifically, occupation number states, used in the quantum interpretation of the os- 
cillations, are found to be much more unstable than coherent ones, typical of the 
semiclassical explanation. 

1 Introduction 

Low-capacitance Josephson-junction devices have recently attracted a wide interest, both 
theoretically and experimentally, particularly in view of the possibility of identifying macro- 
scopic quantum phenomena in their behaviour. In this respect, one of the circuits that have 
gained great attention is the so-called Superconducting Cooper Pair Box (SCB), with an 
increasing number of experiments aimed at supporting a qubit interpretation of its evolution 
(e.g., see PQ-0). 

The SCB is a circuit consisting of two superconducting electrodes linked through a Joseph- 
son junction and coupled, capacitively, to a voltage source. One of the superconducting 
electrodes is assumed to be small enough for the charging energy to play the main role. In 
this situation, there is the possibility of tunnelling electrons one by one through the Joseph- 
son junction, allowing for external control of charge oscillations PQ-[I]. Because of the large 
number of Cooper pairs in the two electrodes, these oscillations have been interpreted as 
genuine (macroscopic) manifestation of quantum coherence. Nevertheless, the experimental 
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data allow for an equivalent explanation of the oscillations in terms of a semiclassical be- 
havior of the system. Two different theoretical scenarios can then be used in modelling the 
observed charge oscillations: one is the so-called quantum phase model [5J E] that essentially 
describes a quantized non-linear harmonic oscillator, while the other is a mean-field model 
that leads to a Gross-Pitaevskii like equation [TUl [11] ; the quantum phase model describes 
the system in terms of occupation number (Fock) states, on the other hand, in the mean-field 
formulations coherent-like states naturally appear. 

More in general, the difficulties for a distinction between classical and quantum behavior 
in Josephson-j unction based devices have been pointed out before in [12]- [Hj. There, it is 
shown that classical nonlinear-oscillator-like models can reproduce some experimental results 
that have been previously attributed to genuine quantum macroscopic behavior. 

In the present work, we limit our consideration to the SCB charge oscillations and ad- 
dress the problem of discriminating between their quantum and classical behavior from an 
open quantum systems perspective, i.e. when the SCB is immersed in a weakly coupled 
external environment. In this case, the SCB dynamics is no longer unitary; instead, it 
is described by a generalized evolution of semigroup type (a so-called quantum dynamical 
semigroup) , that incorporates effects of dissipation and decoherence induced by the environ- 
ment [15]- [18]. More specifically, we will focus our attention on the stability properties of the 
system against the induced environmental noise. A preliminary investigation on this issue 
has been reported in [TJJ], using the so-called singular coupling limit [2U], corresponding to 
a specific environment, with stochastic or white noise correlations. Here instead, we shall 
consider a generic environment and study the SCB stability properties in the weak coupling 
limit |21| . a procedure that implements in a physically consistent way the weak interaction 
of the SCB with the environment. 

The main point of our investigation is that any coupling with an external environment 
needs a microspic description: only in this case one can derive a master equation valid in any 
physical situation |16j . Once a general master equation is obtained, the choice between the 
two possible explanations of the observed SCB charge oscillations depends on the stability 
against noise of the states on which the two different models are based. In this respect, our 
approach differs for instance from the one in [22] , which assumes the validity of the quantum 
phase model and describes the interaction with the environment by means of a spin-boson 
model. 

We found that although both quantum phase and mean-field models predict decoherence, 
they give rise to quite different decay properties: occupation number states (used in the 
former model) turn out to be much less stable than coherent ones (typical of the latter), 
thus confirming the results of [19J. Therefore, the measure of the decay rate of the SCB 
charge oscillations in presence of noise would in principle allow to discriminate between the 
two models. 
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2 SCB oscillations: two approaches 



In a suitable regime and in absence of the environment, the dynamics of a SCB can be 
effectively modelled by a Bose-Hubbard Hamiltonian [231 121] (henceforth, the small electrode 
will be labeled by L, while the other much bigger one by R); in terms of bosonic creation 
and annihilation operators aj, a,, i = L, R in the two electrodes, one can write 

H = E c (a f L a L ) 2 + U L a ] L a L + U R a R a R - K (a j L a R + a L a R ) , (1) 

where the quadratic term Ecia^ai) 2 accounts for Coulomb repulsion in the small island 
(the one in the much larger electrode R can be neglected), L^ajSi, i = L,R, are potential 
contributions, while the last one is the tunneling term. 

Due to the conservation of charge in the SCB, the Hamiltonian ([1]) must be restricted 
to the subspace with a constant total number of particles N = til + tir- Two effective 
descriptions of the dynamics of the SCB can then be obtained, the quantum phase model 
and the mean-field one. 

In the former case, the relevant states are occupation number (Fock) states with a constant 
total number of particles 

|n) = \n L = n,n R = N - n) = - 1 -(a[) n (a R ) N ~ n \v&c) , (2) 

y n!(iv — n)\ 

where |vac) is the vacuum state. Ignoring a constant term, the Hamiltonian in (DO) can be 
rewritten as 

H = E c (n L - n - n g f - K(a[a R + a L a R ) , (3) 

where hi = a) L aL is the number operator in electrode L, while n is the corresponding average 
number; in typical experimental situations, one has n ~ 10 s pQ- [5]. Thus, the difference 
n' = n L — n measures the excess of Cooper pairs (hence of charge) in the same electrode. 
The parameter n g = (U R — Ul)/2E c — n is connected with the potential energy difference 
across the Josephson junction, or equivalently with the average charge in the gate capacitor, 
and it can be controlled externally through the voltage source. 

In the occupation number representation, the Hamiltonian Q takes the form 

N N 

H = E c ^2{n-n-n g ) 2 \n){n\- Ej ^{\n){n + l\ + |n + l)(nf) , (4) 

n=0 n=0 



where we have set K^J [n + l)(iV — n) ~ K^Jn(N — n) = Ej since we are interested in low 
lying states, for which \n — n\ = \n'\ <^ n. Introducing the conjugated operators h' and 0, 
obeying the canonical commutation relations, [0, n'} = i, so that e ±l< ^ decrease (increase) n' 
by 1, the expression in (jlj) is equivalent to the so-called quantum phase model Hamiltonian 

H Q = E c {n - n g ) - Ej cos <p . (5) 
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By adjusting the gate voltage so that n g ~ 1/2, one enters a particular situation (reso- 
nance) where only the occupation number states |n) = \n' + n) for which n' = and n' = 1 
play a role and are strongly coupled by the Josephson junction. In this case, the Hamiltonian 
01]) can then be approximated by a two level Hamiltonian 

Heff = -\[E C {1 - 2n g )a z + Eja x ] , (6) 

where o x and oz are Pauli matrices. The effective two level system will thus display coherent 
oscillation with frequencies given by 



cu q = AE= yjE&l - 2n g f + Ej « Ej , (7) 

where the expression in the middle has been approximated to Ej because we are working 
close to resonance. 

A different description of the system can be given by treating the starting microscopic 
Hamiltonian (jTJ in a mean-field approach. This approximation is justified by the large 
number of Cooper pairs in the two electrodes, all in the same condensed state. This situation 
can be properly described by the product of iV single Cooper pair states 

1 N N 

l*)iV = —7= + Or J |vac) = ^2 C n\ n ) , ( 8 ) 



ra=0 



with 



where ipi and ipn can be interpreted as the condensed wave function in each side of the 
junction (with \iPl\ 2 + \4>r\ 2 = 1)- The dynamics of |^)jv follows the standard Schrodinger 
equation governed by the Hamiltonian (JT]) 

£\*(t)) N = H \*(t)) N . (10) 
at 

The two sides of this equation can be explicitly computed with the help of the following 
relations 

a L \V(t)) N = VN^ L \V(t)) N _ 1} (11a) 

a R Wt)) N = y/N^it))^ (lib) 

n L = N (*(t)\ala L \*(t)) N = N\^ L \\ (11c) 



n R 



N-n L = N (*(t)\tf R a R \*{t)) N = N\ij R \ 2 . (lid) 



By further using a} L <XL ~ (a^ai) = ul = N \ipL{t)\ 2 , justified by the mean-field approxima- 
tion, one finds that the product state |\I/(£))jv is a solution of the evolution equation ([101) 
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if the amplitudes ipi{t) solve the Gross-Pit aevskii equations for the two component order 
parameter (i^Lj^r)'- 



= [U L + NE c \Mt)\ 2 ]Mt) ~ KM*) 
i^ R (t) = U R ip L (t) - Kip L (t) . 



(12) 
(13) 



By setting ipi = ^/rii/ Ne l6i and using the conservation of the total number of particles 
n L + n R — N, the equations above can be written in terms of the Hamiltonian function (up 
to an additive constant) 



H{6,n L ) =E c (n L 



n — n n 



E j cos 9 



(14) 



where the same definitions as before for n g and Ej have been used, together with 6 = 6 l — 9r. 
This Hamiltonian function describes semiclassical charge oscillations with frequency 



<2E c Ej . 



(15) 



Notice that the state \^)n m © behaves like a coherent state; indeed, due to the large 
numbers involved N ^> n ^> 1, we can replace the coefficients C n of \^>)n by a Poisson 
distribution so that in the limit of iV large, one can write: 



N 



n=0 L 



n" 



TV. 



-ind 



n) 



(16) 



a 



n e 



-i0 



which is indeed a coherent state satisfying 

b\a) = a\a), 

where b is the annihilation operator for a fictitious nonlinear oscillator with eigenvectors |n), 

a? L aji\n) = \/N — n b*\n), a^a^Jn) = \/N — n + 1 b\n). 
Hence, the oscillations observed in a SCB might be the result of a semiclassical behaviour de- 



scribed by the coherent state 
at the macroscopic level. 



\a), rather than a manifestation of quantum coherence 



3 SCB with noise: weak coupling limit 

From the previous discussion, it follows that there are two possible scenarios to describe 
charge oscillations in SCB, namely the quantum phase model, which is based on a purely 
quantum description, and the mean-field one which is semiclassical in nature. They both 
are qualitatively consistent with the experimental data: if these could be used to measure 
the oscillation frequency, a direct discrimination between the two approaches would be pos- 
sible; however, the present experimental accuracy does not seem to allow this. Instead, we 
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shall show that the two approaches give different decoherence patterns when the SCB is 
weakly coupled to an environment which acts as a source of noise and dissipation!]] In [TS] 
a specific model of environment has been considered, which required the so called singular 
coupling limit technique; in the following we shall study the effects of a more general source 
of dissipation, using the weak coupling limit approach. 

We describe the weak coupling of the SCB to an external environment by means of the 
total Hamiltonian 

H = H + He + A (ai<& <g> B + a\a 2 <g> £+J , (17) 

where H is the microscopic Bose-Hubbard Hamiltonian (JT]), He is the Hamiltonian of the 
environment, A 1 is a small coupling constant and B is a suitable environment operator. 
We shall assume the environment to be in an equilibrium state p#, with two-point correlation 
functions, (B\t) B) E = Tr e [pe B^(t) B] and similar ones, that decay fast enough (for details, 
see |21j). A heat bath, with Pe — c~^ He , is a typical example of environment fulfilling these 
conditions: specifically, it can be identified with the cloud of non-condensed electrons in the 
two SCB electrodes. 

This very general situation provides the setting for the so-called weak coupling limit [2T] , a 
physically consistent and mathematically precise procedure leading to an evolution equation 
for the SCB density matrix p in Kossakowski-Lindblad form [T5]-|18j: 

^ = -i[H + H^\p]+V[p}. (18) 

The contribution H&> is an environment induced Hamiltonian correction to the starting 
system Hamiltonian Hq, whose explicit expression will not be relevant in the sequel; on the 
other hand, the term D describes the non-Hamiltonian effects of the environment on the 
dynamics of the SCB, which typically result in dissipation and noise. In deriving (|18|) . we 
have assumed to work in the experimentally relevant regime in which the effective Josephson 
coupling Ej is of the same order of magnitude of the dissipative effects, that start to become 
relevant at order A 2 . In this regime, we explicitly find: 

JV 



V W = {*K) (w^n) P W(n) - \{w{n)W\n), pfj + 

+n(u n ) (w(n)pW\n) - \{w\n)W{n), p 



where 
and 



(19) 

2 {w\n)W{n), p)^ , 
W(n) = ^{n + 1)(N - n) \n)(n + 1| , (20) 

/oo poo 
dte'^ {B{t) B^) E , K(u n ) = / dte lt ^ (B\t) B) E , (21) 
-oo J — oo 



1 Within the qubit interpretation, different aspects of decoherence phenomena in SCB behaviour have 
been discussed in [22]. 
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are the Fourier transform of the environment correlations with respect to the frequencies 



u n = E c [2(n - n - n g ) + 1} , (22) 

where n is, as before, the actual number of Cooper pairs in the island lJ§ 

The stability properties of any initial SCB pure state, p = \ip)(ip\, can be studied using the 
master equation (j!8j) . Indeed, first note that the difference between two quantum states can 
be estimated by looking at the part of one of them that is orthogonal to the other. Thus, we 
can quantify the degree of stability of any state p by measuring how fast it starts deviating 
from itself, i.e. by evaluating the initial rate of variation T of the orthogonal contribution 
to itself. Explicitly, 



r= -Tr[(l-|V>>(V>i)p(*)] 



i=0 



dp(t) 
dt 



(23) 



i=0 



where in the last equality we have used (I18p . taking into account that the Hamiltonian 
contribution vanishes. Roughly speaking, a generic state is expected to decay exponentially 
with a rate given by T, so that the bigger T is, the faster it decays, the less stable it is. In the 
specific case under study, we expect to find an appreciable difference between the decay rates 
of occupation number (Fock) and coherent states that might shed light on the suitability of 
the two discussed models. 

Inserting (fl9l) in (|23"1) . in the case of an initial Fock state p = \n)(n\ we get: 



r 



Fock 



{n\T>[\n)(n\]\n} = A 2 (n + 1) (N - n) h{u n ) + n(N - n + 1) k(o;„_i) 



(24) 



while for the coherent like states ([8]) one finds 

Tcohercnt = — <*Jv|X> [\9 N ) (V N \] \9 N ) 
N 

A 2 ]T(n + 1)(N - n)[\C n \ 2 (l - |a + i| 2 ) h(u n ) + \C n+1 \ 2 (l - \C n \ 2 ) k{u u )] . 

(25) 



n=0 



Since N ^> n ^> 1, following arguments similar to those leading to fll6p . the contributing 
terms to the sum can be approximated as |C n+ i| ~ \C n \ ~ 10~ 4 , so that (1 — \C n \) ~ 1; as a 
consequence, one can write 



N 



Tcoherent ~ A 2 J^(n + 1)(JV - n)\C n \ 2 [h(u n ) + K(u n )] . 



(26) 



n=0 



2 There is an additional contribution to T> that arises strictly at resonance, i.e. when n g — 1/2. For 
simplicity, we have omitted it, since it will play no role in the following discussions; indeed, its contribution 
to the decay rates is suppressed by a factor \j\ffi ~ 10 -4 with respect to the dominant one coming from 
(TT9|) (for further details, see [25]). 
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To proceed further, we shall consider a very common instance of environment, that of a 
heat bath having two-point correlation functions of exponentially decaying form, 

(B\t) B) E = (B(t) B^) E = g 2 exp(-\t\/r E ) , 

where t e is the characteristic time scale of the environment and g 2 ~ \(B 2 ) E \ is a constant 
measuring the strength of the bath correlations. Recalling the expression of the frequency 
u n in (122]) . it is convenient to introduce the new integer variable k, by writing n = n + k; at 
resonance, n g = 1/2, one then has u>n+k = 2Eck. The Fourier transforms (12T]) can now be 
explicitly computed, giving 

roc 2 O 2 Tc 

h k = h(u n+k ) = K(u n+k ) = dtg 2 e -W-r B +i^ +k t = y E ; (27) 

1 + (rk) 



3C 



where 

r = — = 2E c r E , (28) 
i~c 

with re = (2Ec)~ l being the characteristic time of oscillations due to Coulomb interaction^ 
With these results, using the central limit theorem to approximate ICV^I 2 by a Gaussian 
distribution^ r coherent in f[2"6T) can be cast in the following form: 

N-n „ 9 

n x — \ 2 o Tp e 2n 

Inherent « A 2 V (fi + k + 1)(N - 71 - k) ~ \ ■= , (29) 

k t^n l + {rk) V2irn 

that, in turn, due to the large numbers involved, is well approximated by 

i r°° e~ y2/2 

Tcoherent « A 2 (2 g 2 t e ) n(N - n) f(Vfir) , f(z) = -= / Ay- . (30) 

V27T 7-oo l + (^y) 

One can similarly evaluate the expression of the decay rate in the case of occupation number 
states, obtaining from (pMl) : 

r Fock ^ \ 2 {2g 2 T E )n{N -n) . (31) 

The integrand in f(z) above is a product of two decaying functions with different character- 
istic scales, namely a Gaussian with standard deviation equal to one and a Lorentzian with 
width 1/z. Thus, in the large z limit the integral is dominated by the Lorentzian. Hence, 
for y/Er ^> 1, we get 

' '" k oc \ffir , (32) 



coherent f ( \pWT 



3 In the experiment reported in [HE], the charging energy is around Ec ~ 500 fieV, so that tc ~ 10 11 s. 

Substituting (Tile]) . (fTTd]l in the definition ©, one finds |C„| 2 = n! (#l ra)! (f )" (l ~ f ) N ~ n , since 
(ni) « n. Using the Stirling formula for the various factorials, one further gets |C R+ fc| 2 w '^^■\J^ N-n) g(k) 
where, in the limit N n, the function <?(fc) can be very well approximated by a Gaussian distribution (for 
further details, see [25"]). 
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showing that the decay rate in the case of occupation number states results much larger than 
the corresponding one for coherent states. 

However, to be experimentally relevant, this result needs to be interpreted within the 
actual conditions of a typical setup; in particular, one needs to estimate the magnitude of 
these decay rates and compare them with the smallest, characteristic energy scale of the 
device under study, i.e. the tunneling energy Ej. 

Indeed, as already mentioned, in order for the noise effects to be observable, the strength 
of the interaction describing the coupling of the SCB with the environment should be of the 
same order of magnitude of the tunneling term. Recalling the form of the corresponding 
Hamiltonian pieces (pQ) and (TTTj) . this condition roughly means or equivalently 

Xg^n(N — n) m Ej. Then, from the expression fl3T|) . one immediately gets the estimate: 

TFock _ Ej . . 

-ET- ~-FT r ■ (33) 
tjj hj C 

The ratio Ej/Eq is fixed in any experimental situation. For instance, in the setup described 
in [H [2] , one has 

Ej « 50 /xeV » 10 10 s _1 , 
E c « 500 //eV » 10 11 s _1 , 

so that Ej/Ec ~ 1/10. As a consequence, in the relevant regime \ffir 3> 1, the decay rate 
Tpock might be close to Ej and therefore the effects of the environment visible. In particular, 
by choosing an environment for which r is close to onej^l one obtains T-p OQ k ~ Ej/10, or 
equivalently a decay time of order r Fock = l/r Foc k ~ 10/ Ej « 10~ 9 s. 
On the other hand, for the coherent states one gets 

1 Ej (34) 



Ej y/nEc 

which is independent of r and, in the case of the setup in [TJ [2], yields a decay time of 
order r coherent = l/r cohcrcnt ~ 10~ 5 s. This means that if the charge oscillations are a 
manifestation of macroscopic quantum effects, their damping would become evident after 
just few oscillation periods, while if they are semiclassical in nature, they will persist for 
very long times. In the specific case of the experiment described in [HE], by doubling or 
tripling the observation time in presence of a suitably engineered environment, damping of 



5 This estimate comes from imposing (Hj) ~ (Hi), where Hj is the tunneling term in (jTJ) and Hi is 
the interaction term in (fT7|) . Indeed, (Hj) ~ K^lo^), while (XHi) ~ X(B)(a,La r R ) ~ Xg^LO^), since 

6 This choice is compatible with the weak coupling limit. As mentioned before, the application of this 
procedure requires that the decoherence time scale of the slow subsystem dynamics r = 1/r be much greater 
of both the characteristic decay time te of correlations in the environment and the intrinsic time scale of the 
free system dynamics tc = l/2Ec- One checks [25] that with the choice r ~ 1 both conditions are satisfied. 
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the charge oscillations must be observed if the oscillations are indeed a manifestation of 
quantum coherence. 

As a final remark, note that the ratio u c /u q between the charge oscillation frequencies 
in (JT5j) and (jTJ) behaves as (Ec/Ej) 1 / 2 ; thus, by modifying the experimental conditions so 
that to increase the ratio Ec/Ej, one can eventually distinguish between the two types of 
explanation of the observed charge oscillations through a frequency measure. Nevertheless, 
the alternative approach here proposed of adding noise to the device might be, in principle, 
more efficient in discriminating between the two models. Indeed, the damping time for Fock 
states increases linearly with Ec/Ej, while the ratio of the frequencies above only with the 
square root. 

4 Conclusions 

Reviewing the phenomenology of charge oscillations in a superconducting Cooper pair box, 
we have seen that they might be modelled in two different ways, either as a manifestation of 
macroscopic quantum coherence (leading to the widely accepted qubit interpretation) or as 
the result of a semiclassical mean- field approach (through a Gross-Pitaevskii like equation). 
However, the response of the SCB to external noise, generated by a weakly coupled environ- 
ment, is found to be very different in the two models. Both approaches predict damping of 
the oscillations, but the decay rate in the case of qubit (Fock) states differs by a factor \ffir 
from that of the mean-field (coherent like) ones. In the physically relevant regime \/Er ^> 1, 
one then expects the Fock states to decay much faster than the coherent ones. 

This result might provide a way to distinguish experimentally between the two possible 
interpretations of the observed SCB charge oscillations. The idea is to couple an SCB with 
an externally controlled environment, satisfying the conditions of the weak coupling limit 
approximation. As soon as the interaction with the environment is switched on and therefore 
noise is injected into the nanodevice, damping of the charge oscillations should become visible 
if these are a manifestation of macroscopic quantum coherence, as described by the qubit 
model; on the other hand, if they survive for long times, this would be an indication of their 
semiclassical origin, as described by the mean-field approach. 

One can easily evaluate the visibility of the damping effects in the specific setup described 
in [11 [2]. For the physically relevant case of a heat bath with r ~ 1, we have found that 
the decay time for occupation number, Fock states can be estimated in about 10~ 9 s. As 
a consequence, if due to macroscopic quantum coherence, the decoherence effects must be 
visible after about a few oscillation periods. In the case of coherent states, we have instead 
obtained a much longer decay time of about 10 -5 s. Therefore, the possibility to experi- 
mentally distinguish between the two types of interpretation with already existing setups 
appears quite realistic. 
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